Rate of Period Change as a Diagnostic of Cepheid Properties 
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ABSTRACT 

Rate of period change P for a Cepheid is shown to be a parameter that is capable of indicating 
the instability strip crossing mode for individual objects, and, in conjunction with light amplitude, 
likely location within the instability strip. Observed rates of period change in over 200 Milky Way 
Cepheids are demonstrated to be in general agreement with predictions from stellar evolutionary 
models, although the sample also displays features that are inconsistent with some published 
models and indicative of the importance of additional factors not fully incorporated in models to 
date. 



Subject headings: stars: Cepheids — stars: evolution 



1. Introduction 

Cepheids represent a brief phase in the post- 
main-sequence evolution of stars t hat originall y 
had masses in excess of - 3±M riTurnerlll996h . 
As evolved objects they populate the Cepheid in- 
stability strip in the HR diagram according to the 
manner in which they generate energy, depending 
upon strip crossing. Intermediate mass stars can 
become unstable to radial pulsation during shell 
hydrogen burning (first crossing), twice during 
core helium burning (second and third crossi ngs), 
and according to some evolutionary mo d els (llben 
1965[ iBecker et all Il977t iBeckerl Il985t IXu fc Li 
2004), twice during shell helium burning (fourth 
and fifth crossings). A common feature of many 



recent evolutionary models for intermediate-mass 
stars, which incorporate the new opacity tables 
/e.g. 
l20Qi 



0; 



wnicn incorporate trie new opacity tables 
Mevnet fc Maederl |2000L l2002t iBono et al 



Salasnich et al.l 



2000), is that they permit 



only three strip crossings for such stars, since core 
oxygen ignition occurs prior to a separate shell 
helium-burning phase. 

In all cases the evolution of stars through the 
Cepheid instability strip should be associated with 
gradual changes in overall dimensions, and hence 
periods of pulsation: period increases for evolu- 
tion towards the cool edge of the instability strip 
as the stellar radius grows, and period decreases 
for evolution towards the hot edge as the stellar 
radius decreases. The observed parabolic trends 
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in Cepheid 0-C diagrams (plots of the differences 
between Observed times of light maximum and 
those Computed from a linear ephemeris) have 
been recognized for the past half century as ev- 
idence for the evol ution of such stars through 
the instability strip dParanego 19581 : IStruvej 1959; 
lErleksova fe Irkaevl Il982 ). As noted by IStruvd 
(|l959h . "It appears that studies of period change 
are by far the most sensitive test available to the 
astronomer for detecting minute alterations in the 
physical characteristics of a star." 

Observations of period changes in Cepheids 
have been matched with some confidence to evo- 
lutionary models of massive stars in various cross- 
ings of the instability strip (e .g., iTurner 1998; 
ITurner fc Berdnikovl l2001. 2004) in order to iden- 
tify the direction of strip crossing for individual 
variables. When used for such purposes, the study 
of Cepheid period changes becomes an important 
tool for the characterization of individual members 
of the class. 

In principle it should also be possible to use 
rate of period change for individual Cepheids to es- 
tablish likely location within the instability strip. 
Because strip crossings for individual Cepheids oc- 
cur at different rates and at different luminosi- 
ties for specific stellar masses, the observed rates 
of period change must be closely related to strip 
crossing mode and location within the instabil- 
ity strip. Potential constraints are imposed by 
variations in chemical composition and pulsation 
mode, e.g., fundamental mode, first overt one, etc. 
(|Berdnikov et al.lll997t ITurner et alll999h , as well 
as by our limited ability to establish small rates of 
period change for O-C dat a containing siz eable 
observational uncertainties ( Szabadosl 1983). In 
this paper we demonstrate the link in more detail. 

2. Basis of the Relationship 

The link between rate of period change in 
Cepheids and location within the instability strip 
is illustrated with the aid of Fig. 1. The dia- 
gram is a theoretical HR diagram that depicts 
the location of the Cepheid instability strip ac- 
cording to the paramet ers derived for Milky Way 
Cepheids (jTurnerl 120011 ). along with Geneva evo- 
lutionary tracks for stars of 4, 5, 7, and lOM© 
at Z = 0.008 from iLeieune fc Schaererj (|200l[ ). 



various portions of the instability strip. Accord- 
ing to the well established Cepheid period-radius 
relation, they should represent lines of constant 
pulsation period for individual Cepheids. 

^From an examination of Fig. 1 it is clear that, 
if one considers only Cepheids of a specific period 
and in a common crossing of the instability strip, 
those on the hot edge of the strip must be ~ 20% 
more massive than those on the cool edge of the 
strip. Since rate of evolution increases in propor- 
tion to the mass of a star, Cepheids lying on the 
hot edge of the strip are evolving faster, and hence 
changing their pulsation periods at a more rapid 
rate, than Cepheids of identical period lying on 
the cool edge of the strip. Rate of period change 
therefore relates directly to location within the in- 
stability strip for individual Cepheids. Differences 
in strip crossing modes are only a minor concern. 
Cepheids with increasing periods must be in the 
first, third, or fifth crossing of the strip, whereas 
Cepheids with decreasing periods must be in the 
second or fourth crossing of the strip. 

A minor complication arises from restrictions 
on our ability to identify period changes in 
Cepheids tied solely to stellar evolution. Some 
Cepheids exhibit erratic period changes that ap- 
pear to originate from r andom fluctuations in pul- 
satio n period. SZ Tau (Berdnikov fc Pastukhoval 
1995h. S Vul (lBerdnikovl ll994h. and V1496 Aal 
(jBerdnikov et all |2004) ~ are excellent examples, 
although in the first two cases it is possible to 
identify the underlying evolutionary modifications 
to pulsation period. 

A study by Berdnikov fc Ignatova (2000) may 
give the impression that stellar evolution has only 
a minor effect on Cepheid O-C diagrams, since it 
notes that parabolic trends were detected in only 
67 of 230 Cepheids surveyed. That number is mis- 
leading, howe ver, given that a previous survey by 



Lines of constant stellar radius are shown crossing 



Turner! (|1998l ) had found parabolic trends in 137 
Cepheids from a much smaller sample. It was 
actually intended to indicate the poor temporal 
coverage and lack of extensive O-C data available 
for many well-studied Galactic Cepheids, a situ- 
ation that has been remedied in recent years by 
our ongoing program to obtain archival data on 
Cepheid brightness variations using the Harvard 
College Observatory Photographic Plate Collec- 
tion. At present the parabolic trends in O-C dia- 
grams typical of stellar evolution are found to be 
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extrem ely common. A survey bv lGlushkova et al 
( 20051 ) cites a typical frequency of ~ 80% in both 
cluster and field Cepheids, for example, although 
their "anomalous" objects include Cepheids like 
SV Vul in which the evolutionary trend is quite 
distinct (jTurner fc Berdnikovll2004r ). A more real- 
istic frequency for Milky Way Cepheids displaying 
evolutionary trends is in excess of ~ 90%. For 
many of the remaining objects, the evolutionary 
trends may be more obvious in longer time base- 
lines of light curve coverage. 

As also pointed out by Fernid ( 1990h and by 
Berdnikov fc Turnerl (|2004h . the O-C trends in- 
dicative of evolution in Cepheids need not be 
strictly parabolic. If the rate at which a mas- 
sive star is evolving through the instability strip 
is not constant with time, the O-C data for the 
associated Cepheid variable may be better de- 
scribed by a third or fourth orde r polynomial. 
The Cepheids Y Oph (lFernid ll990h and WZ Car 
( Berdnikov fc Turner 2004 ) are two objects (of 
several hundred) where that appears to be the 
case. Such complications may affect the derived 
rates of period change, but in most cases only by 
small amounts. In the large majority of studies of 
Cepheid period changes, the derived rate of period 
change reflects the evol ution of the sta r through 
the instability strip (see Szabadosl 1983f ) . 



3. Stellar Evolution Predictions 

Most computational evolutionary models for 
evolved stars are used for constructing evolution- 
ary tracks rather than testing for pulsation insta- 
bility. But, as noted by citetpa58, it is possible 
to use the basic information they provide on grad- 
ual changes in luminosity and effective tempera- 
ture to predict expected rates of period change for 
Cepheids of different period. A starting point is 
the well known period-density relation: 



PM 2 



= Q 



where P is the pulsation period, p is the density, 
M is the stellar mass, R is the stellar radius, and 
Q, the pulsatio n constant, has a small p eriod de- 
pendence (e.g., Kraft 1961 : Fernid 1967 ) that we 
assume here varies as P» bas e d upo n an empirical 
analysis bv iTurner fc Burkd ([2002). Differentia- 
tion of the period-density relation, in conjunction 



with the standard equation for stellar luminosity, 
therefore leads to the following result: 



P 
P 



6L 24 T 
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The desired quantity, the rate of period change P, 
is obtained from tabulated differences in stellar lu- 
minosity and effective temperature as a function of 
age as a model star evolves through the instability 
strip. 

For the present study we calculated values 
of P from the above relationship using com- 
putational stellar evolutionary models from a 
variety of availab l e pub lis hed 



iMaeder fc Mevnetl d 19881). lAlibert et al 



sources, namel y 



Claret 



1999), 



2004). 



iLeieune fc Schaererl and 
The published data were used to compute dif- 
ferent parameters, depending upon the availabil 
of the 



i ty of the necessary information. lAlibert et al 
(1 19991) cite parameters for stars of different mass 
reaching the hot and cool edges of the instabil- 
ity strip, so their data yield information only 
about rates of period change nea r the center o f 
the strip. In other cases, such as IClaret (2004). 
there is sufficient time resolution in the output 
parameters to track changes in pulsation pe- 
riod across individual instabil ity strip crossings 
For the remaining sources (jMaeder fc Mevnet 



1988 : Lei eune fc Schaererl 20011 ). including IClaret 
(120041) . we calculated rates of period change for 
the intersection of the evolutionary tracks with 
the observationally delineated boundaries of the 
i nstab ility strip defined empirically by ITurner 
(2001), which are close to t hose predicted by m od- 
els of pulsation instability ( Alibert et al.lll999h . as 
well as for points lying within the strip bound- 
aries. Pulsation periods were established usin g 
the period-radius relation ( Turner fc Burkdl2002j ). 
The present results differ from those obtained ear- 
lier (jTurner fc Berdnikov|[200lll2003l) in being tied 
to a larger variety of models with a greater range of 
metallicity, and by the inclusion of a weak period 
dependence for Q in the period-density relation. 

The computed results on rates of period change 
are plotted in Fig. 2 for all of the accessible mod- 
els. Different symbols denote the di fferent sources. 



Values calculated from the models of Alibert et al 



(1999) are plotted using filled circles, while others 
are plotted using open circles. Plus signs indi- 
cate results calculated for stars evolving through 
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the hot and cool edges of the instability strip, with 
the rate of period change in general being larger on 
the hot edge of the instability strip, i.e., for more 
massive stars. Large symbols denote stars of solar 
metallicity, Z — 0.02, intermediate-sized symbols 
denote stars with metallicities of Z = 0.01 and 
Z = 0.008, and small symbols denote stars of very 
low metallicity, Z = 0.001 and Z = 0.004. Lines 
have been drawn to enclose those regions within 
which the results for different crossing modes ap- 
pear to cluster. Sequences of points indicate mod- 
els for which the time resolution was fine enough 
to calculate rate of period change over the entire 
crossing of the instability strip. 

The distribution of data points in Fig. 2 sug- 
gests a variety of different conclusions regarding 
the models. First, the different models for the 
rapid first crossing of the instability strip are in 
very good agreement, and display very little varia- 
tion with metallicity. The first crossing of the strip 
is a rapid transition for all stars, regardless of in- 
dividual differences in rotation rate, etc., and that 
is evident from the models. Evidently the com- 
putational codes used for calculating the phases 
of shell hydrogen burning in stars, while perhaps 
differing in detail from one source to another, gen- 
erate nearly identical results, the small variation 
in rate of period change at specific pulsation pe- 
riod arising from the finite width of the instability 
strip and the fact that more massive stars cross 
the strip at a greater luminosity and at a faster 
rate than less massive stars. For stars in the first 
crossing of the strip, high rate of period increase 
at specific pulsation period corresponds to stars 
on the hot edge of the strip, low rate of period 
increase to stars on the cool edge of the strip. 

Negative period changes arise during the second 
crossing of the instability strip, which occurs dur- 
ing the blue loop phase of stellar evolution follow- 
ing the onset of core helium burning. The extent of 
the blue loop can d epend upon a variety of factor s 
(see, for example, iBeckerl 11981 IXu fc Lil [2004ft . 
such as metallicity, the treatment of core over- 
shooting, and the distribution of CNO elements 
throughout the star. All factors affect how far 
a star enters the instability strip during core he- 
lium burning, and presumably affects how rapidly 
it evolves within the strip. Given the potentially 
large differences in initial conditions for such stars 
as main-sequence objects, for example, large vari- 



ations in initial rotation rate, one might expect 
real stars to display large variations in how far 
they penetrate the Cepheid instability stripm as 
core helium burning objects. Somewhat unexpect- 
edly, there are also very large variations among the 
models stars as well. 

Evidently, metallicity plays only a minor role 
in governing the rate at which stars traverse the 
instability strip. There is as much dependence on 
the specifics of the stellar evolution ary code used. 
The models of Alibert et al. (|l999ft . for example, 



generate faster rates of period decrease than do 
other models, despite the use of common opacity 
tables. Models from individual sources are at least 
internally consistent in their predictions for stars 
of different masses and for stars in all portions 
of the second strip crossing. The rates of period 
decrease during individual strip crossings are also 
very similar to the variations predicted on the ba- 
sis of mass differences, i.e, predicted variations in 
rate of period decrease at a specific pulsation pe- 
riod are generally small, except for long period 
Cepheids. 

The third crossing of the instability strip oc- 
curs during the late stages of core helium burn- 
ing, and gives rise to period increases, for which 
the predicted rates are depicted in the top portion 
of Fig. 2 along with those for the first crossing. 
Most of the comments regarding the second cross- 
ing of the strip apply equally to the third crossing. 
Again, metallicity seems to play a less important 
role in the predicted rates of period increase than 
di fferences in the evolu tionary code. The models 
of Alibert et al.1 ( 1999ft predict faster rates of pe- 
riod change (period increases in this case) than 
do other models, although with less consistency 
for stars of different mass. The rates of period 
increase during individual strip crossings are also 
similar to the variations predicted on the basis of 
mass differences, and predicted variations in the 
rate of period increase at a specific pulsation pe- 
riod are generally small. 

A well-known problem arises for low-mass stars 
in the second and third crossings of the instability 
strip, since the blue loop phases of evolutionary 
models for stars of solar metallicity, Z = 0.02, 
do not enter the strip for M < 4.75M (see 
I Alibert et al. 1999ft . Model stars of lower metal- 
licity can traverse the strip at smaller masses, 
but often only on the cool edge. By inference, 
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most classical Cepheids of near-solar metallicity 
should ha ve pulsation p eriods in excess of ~ 3^ 
days fe.g.. lTurnerlll996j ). consistent with the ob- 
servational sample. Nearby Milky Way Cepheids 
have abundances close to the solar values (e.g., 
Andrievskv et al] l2002al lbh . and only a few have 
periods of less than 3^ days. Many may be over- 
tone pulsators. 

The observational picture is illustrated in 
Fig. 3, which presents available data on pe- 
riod changes for over 2 00 Cepheids, as obtained 
from t he literature ( Berdnikov fc Pastukhova 
1994allbl.[l995l:lBerdnikov et al]ll997tlTurnerll998l: 



Berdnikov fc Ignatoval200d : IBerdnikov et alf2 003) 
and ongoing research studies by the authors (e., 
Berdnikov fc Turner|[200l iBerdnikov et al.llm 



The relationships plotted in Fig. 3 depict the re- 
gions within which the model calculations appear 
to cluster. 

It has been pointed out previ ously (e.g., Szabadosl 

19831 : lFerniJl984 lTurnerlll998l ) that the observed 
rates of period change in Cepheids are generally a 
good match to predictions from stellar evolution- 
ary models. The data of Fig. 3 provide further 
confirmation of that conclusion. Moreover, three 
further conclusions can be reached. First, once 
consideration is taken of the expected changes 
arising from evolution through the instability 
strip, the observed period changes in Cepheids are 
unlikely to contain any sizable component arising 
from another source. There are only a few excep- 
tions to such a conclusion, and they are rather 
unus ual objects like V1496 Aql (IBerdnikov et al. 



20041 ). which exhibits period changes dominated 
by random fluctuations in pulsation period. 

Second, the observed period changes in Cepheids 
deviate in small but important ways from what is 
expected according to predictions based upon spe- 
cific stellar e v olutio nary models. The models of 
Alibert et al. ( 1999h . for example, predict faster 
second and third crossings of the strip than those 
observed, and at much different rates. In contrast, 
the observed period changes in Cepheids are very 
similar for objects likely to be in the second and 
third crossings. The models of IClaretl (|2004T ) are 
most consistent with observations in that regard, 
but it is necessary to have a more complete mass 
grid of models constructed in the same manner to 
make a more detailed comparison. 

Third, the range in observed rates of period 



change for most Cepheids is smaller than that 
resulting from a comparison of the results from 
different evolutionary models. That is somewhat 
surprising, given our previous discussion about 
potentially wide variations in initial conditions 
for Cepheid predecessors. Evidently real stars 
are similar enough in their internal characteristics 
that they evolve at fairly similar rates through the 
Cepheid instability strip. 

The proportions of Cepheids in different cross- 
ing modes and in different period ranges in Fig. 
3 are also reasonably consistent with evolutionary 
expectations. For example, stars in the first cross- 
ing of the instability strip during shell hydrogen 
burning are evolving about two orders of magni- 
tude faster than stars in second and third cross- 
ings, so their relative numbers should be small. 
The two Cepheids in Fig. 3 undergoing large rates 
of period increase and falling in the predicted re- 
gion for first crossers are Polaris (a UMi) and DX 
Gem. We assume that both arc first crossers, as 



was also argued for Polaris bv lTurner et al.l (|2005l ). 
Moreover, the observed rate of period change for 
Polaris is now seen to be exactly what stellar evo- 
lutionary models predict for a star lying on the 
cool edge of the instability strip for first crossers. 

The proportion of Cepheids with detectable 
parabolic trends in their O-C data also increases 
noticeably towards short pulsation periods, which 
is again consistent with the evolutionary expec- 
tation that the most abundant pulsators must be 
those evolving most slowly through the instability 
strip. There is a curious anomaly in the distribu- 
tion of short period Cepheids, where essentially no 
variables are found to have rates of period change 
as predicted for stars in second and third crossings 
of the strip at P < 3.5 days (log P < 0.55). Such 
stars have pro genitor masses of less than ~ 4 Mq 
( Turner 1996), where stellar evolutionary models 
for solar metallicity stars predict that the evolu- 
tionary tracks for core helium burning stars should 
no longer enter the strip. The short period cutoff 
in the observational sample is therefore consistent 
with expectations from stellar evolutionary mod- 
els. But the existence of stars of P < 3.5 days 
with rates of period change roughly an order of 
magnitude faster than predicted for stars in sec- 
ond and third crossings of the instability strip is 
not. The uncertainties in the observed rates of pe- 
riod change in the anomalous objects are generally 
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much too small to resolve the anomaly by invoking 
systematic errors in the values of P. 

An additional factor that can be important 
for short period Cepheids is overtone pulsation. 
The Cepheids in the observational sample have 
all been assumed to be fundamental mode pul- 
sators, and require a displacement of +0.15 in log 
P to establish their proper locations in Fig. 3 
if they are overtone pulsators. Yet the applica- 
tion of such corrections to all of the anomalous 
objects does not affect their distribution signifi- 
cantly; most still fall outside the region of P-space 
predicted for stars in the second and third cross- 
ing of the instability strip. Current stellar evolu- 
tionary models are therefore unable to explain the 
existence of such stars, which suggests that the 
manner of treating the details of stellar evolution 
during blue loo p stages is very important (see also 
Xu & Li 2004). That is one area where improve- 
ments to the observational sample on Cepheid pe- 
riod changes can play an important role in testing 
the results from stellar evolutionary models. 

4. P as a Fundamental Parameter 

In Fig. 3 the dispersion in the rates of pe- 
riod change P observed in long period Cepheids 
(P > 10 d ) is smaller than what is observed for 
the calculated dispersion in that parameter among 
different stellar evolutionary models. One might 
expect P to correlate closely with location in the 
instability strip for Cepheids in all strip crossings, 
according to the results of Fig. 1. It is informative 
to examine the observational data more closely to 
determine if that is the case. 

As a first step, we note that the observed rates 
of Cepheid period change plotted in Fig. 3 fall 
mainly within specific bands delineated by linear 
margins of slope 3.0 separated by an order of mag- 
nitude range in P. Fig. 4 is a separate plot of the 
data that displays such empirically defined mar- 
gins. All but two of the long period Cepheids with 
increasing periods fall within the lower set of mar- 
gins, as do the majority of short period Cepheids 
with increasing periods. Cepheids with decreas- 
ing periods display a greater dispersion in P that 
may be intrinsic, or may be caused by larger un- 
certainties in P for the stars, particularly those 
with small rates of period change. 

The anomaly for Cepheids with P < 3.5 days 



(log P < 0.55) is again apparent in Fig. 4. All 
Cepheids of shorter period display faster rates of 
period change than is typical of variables popu- 
lating the lower band, and there are a number of 
stars of longer period also falling in this region of 
rapid period change. Presumably those stars rep- 
resent Cepheids in fourth and fifth crossings of the 
instability strip, with faster associated rates of pe- 
riod change. Multiple crossings of the instability 
strip appear to be possible for stars in late core 
helium burning stages, depending upon the CNO 
abun dances in the hy drogen burning shells of such 
stars dXu fc Lill200"l . 

The finite range in stellar surface tempera- 
ture for stars populating the instability strip at 
constant pulsation period implies distinct differ- 
ences in pulsation efficiency that should coincide 
with marked differences in pulsation amplitude for 
Cepheids of similar period. On the hot edge of the 
strip the ionization zone is just beginning to reach 
depths where the piston mechanism for pulsation 
becomes efficient, so light amplitudes should be 
small but increasing with decreasing surface tem- 
peratures. On the cool edge the lower surface 
temperatures are associated with increased con- 
vective energy transport in the star's outer layers 
(jDeupreel fl980). so pulsation amplitudes should 
also be small. 

The first study of Cepheid am plitude s as a 
function of position in the strip by iKraftl i|l963l ) 
was consistent with that picture, although small 
amplitude Cepheids were found only on the hot 
edge of the strip. All sub sequent amplitud e 
maps of the instability strip b y Hofmcistcr (1967), 



(Il974h. IPel fc Lubl [l978^ lTurneil (l200ll). and 



ISandaee fc Tammannl (fl97lh . Pavne-Gaposchkml 



Sandage et alj (|2004l ) have produced similar re- 



sults, namely a sharp rise to maximum amplitude 
on the hot edge of the strip followed by a more 
gradual decline towards the cool edge. 

Cepheid amplitudes display a period depen- 
dence as well as a dependence upon location within 
the strip, a natural consequence of an effect tied 
to surface gravity as well as pulsation efficiency. 
In order to eliminate that factor in characteriz- 
ing Cepheid period changes, we have normalized 
the resulting values of blue light amplitude and 
P as follows: (i) blue amplitudes AB were stan- 
dardized through the ratio AP/AP(max), where 
AP(max) is the maximum value of AB for the 
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star's pulsation period, and (ii) P was adjusted to 
the equivalent value for a Cepheid with a pulsa- 
tion period of 10 d using the empirically-obtained 
slope plotted in Fig. 4. 

Fig. 5 plots such data for Cepheids with 12 d < 
P < 40 d and increasing pulsation periods (P ~ 
20 d ). The upper part of the diagram plots the 
individual data, while the middle part of the dia- 
gram plots running means for the data. The lower 
part of the diagram is an alternate interpretation 
of the same data, as described below. Similar plots 
are given in Fig. 6 for Cepheids with 4 d < P < 8 d 
and increasing pulsation periods (P ~ 6 d ), and in 
Fig. 7 for Cepheids with 4 d < P < 8 d and de- 
creasing pulsation periods (P ~ 6 d ). Recall that 
large values of P should correspond to the hot side 
of the instability strip, and small values to the cool 
side. 

The data for 20 d Cepheids (top portion of Fig. 
5) display a tendency for large amplitude Cepheids 
to have rates of period increase typical of stars ly- 
ing near the center of the instability strip, with 
smaller amplitude Cepheids falling towards the 
hot and cool edges (larger and smaller values of 
P, respectively). The trend is more obvious when 
one plots running five-point means of the same 
data, as in the middle section of Fig. 5. There are 
two long period Cepheids with anomalously large 
values of P, SZ Cas and AQ Pup, which are con- 
ceivably fifth crossers. If they are omitted from the 
running means and averages over smaller samples 
are included at the extremes of P, one obtains the 
results in the lower portion of Fig. 5, which are 
typical of independent cross-sectional amplitude 
maps of the instability strip. The scatter in P val- 
ues evident in the top portion of Fig. 5 is intrinsic 
to the stars, and is not the result of large uncer- 
tainties in the calculated values. Presumably there 
are intrinsic physical differences from one Cepheid 
to another that account for the scatter, as noted 
earlier. Differences in initial rotation velocity for 
the progenitor main-sequence stars might be the 
sole factor, given that they would generate suffi- 
ciently large variations in the abundances of the 
CNO elements throughout th e star to affect the 
extent of the blue loop stages (|Xu fc Lil 120041 ) . 

The data for 6 d Cepheids with period increases 
(top portion of Fig. 6) are more complicated. It 
appears that the sample consists of two overlap- 
ping groups of objects, a feature that also appears 



in the running five-point means displayed in the 
middle section of Fig. 6. We assume that each 
group consists of Cepheids displaying an order of 
magnitude (factor of 10) variation in P values from 
the hot to cool edges of the instability strip, as dis- 
played by the long period Cepheids in Fig. 5, and 
use the results presented in the lower portion of 
Fig. 5 as a template for the likely variations in rel- 
ative amplitude with P for short period Cepheids. 
When the two groups in Fig. 6 are separated in 
such fashion and averages over smaller samples are 
included at extreme values of P for each group, one 
obtains the results in the lower portion of Fig. 6. 
The simplest explanation for the existence of two 
groups among the short period Cepheids is the ex- 
istence of higher strip crossings among the stars, 
namely fifth crossings for Cepheids undergoing pe- 
riod increases. 

Similar results apply to the data for 6 d Cepheids 
with period decreases (Fig. 7), when analyzed in 
similar fashion, despite the smaller sample size. 
The top portion of Fig. 7 displays excessive scat- 
ter, much like that in the top portion of Fig. 6, 
with only marginal improvement through running 
five-point means (middle portion of Fig. 7). Re- 
stricting the data sets as above produces the re- 
sults depicted in the lower portion of Fig. 7, which 
suggests an overlap between Cepheids in second 
and fourth crossings of the instability strip. As 
noted for long-period Cepheids, there is an intrin- 
sic scatter in P values that is not the result of large 
uncertainties in the calculated values. Such scatter 
makes it difficult to use rate of period change for 
individual Cepheids to identify their exact location 
within the instability strip, although approximate 
placements are possible in most cases. 

The conclusions reached here, while admittedly 
speculative, provide a reasonable explanation for 
the characteristic behavior of pulsation amplitude 
and rate of period change for individual Cepheids 
at specific pulsation period. As noted earlier, 
there may be further complications arising from 
the presence of overtone Cepheids in the sample, 
but their numbers should be relatively small in 
the present case, except at short periods, and they 
would not alter the observed distribution of data 
points significantly. A more comprehensive study 
including recognized overtone Cepheids should be 
possible once O-C studies have been completed 
for such variables. 
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5. Discussion 

Our intent here is to demonstrate that rate of 
period change for a Ccpheid is a useful parame- 
ter that permits one to characterize the variable 
in terms of specific evolutionary state. Informa- 
tion on P for a Cepheid, in conjunction with its 
known pulsation period and light amplitude, can 
be used to identify the strip crossing mode for the 
object as well as its likely location within the strip, 
the latter independent of its observed color and 
reddening. The parameter P may even be use- 
ful for establishing if a Ccpheid is a fundamental 
mode pulsator or an overtone pulsator, although 
we leave that as a future exercise. 

If we interpret the results of Figs. 5-7 as 
a generic indicator of how pulsation amplitude 
varies across the instability strip, then the width of 
the strip in P at constant period, which amounts 
to ~ 1.2 in log P, must encompass a range of ~ 16 
in P. Of that, an intrinsic dispersion in log P 
values amounting to perhaps 0.4-0.5, a factor of 
<~ 3, presumably arises from actual internal dif- 
ferences in the Cephcids resulting from different 
histories for their progenitor stars. Specific stellar 
evolutionary models presented in Fig. 2 predict a 
smaller variation in log P than what is observed, 
which may reflect the simplicity of the models. In 
that regard, observed rates of period change in 
Cepheids can play an important role as a check on 
how closely stellar evolutionary models match real 
stars. Until now Cepheid period changes have not 
been used for that purpose. 
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Fig. 1. — The theoretical HR diagram illustrat- 
ing post-main-sequence evol utionary tracks for 
stars of 4, 5, 7, and 10 M Q (jLeieune fc Schaerer 
200 ll ). Included is the observational location of 



the Cepheid instability strip (dotted lines, from 



Turner 2001) and lines of constant stellar radius. 
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log P(d) 

Fig. 2. — Predicted rates of period change for stars 
crossing the Cepheid instability strip as tied to 
published stellar evolutionary models. The mean- 
ing of the different symbols is explained in the 
text. Lines denote regions within which the pre- 
dictions from different stellar evolutionary models 
appear to cluster. The different crossings of the 
instability strip are identified. 
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Fig. 3. — Observed rates of period change, 
along with their calculated uncertainties, for well- 
studied Cepheids possessing many years of O-C 
data. Dotted lines are the relations depicted in 
Fig. 3, and the different strip crossings are identi- 
fied. 
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Fig. 4. — The data of Fig. 3 plotted along with 
suggested empirical delineations of the regions cor- 
responding to Cepheids in different crossings of the 
instability strip, as identified. 
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Fig. 5. — Normalized blue amplitudes of Cepheids 
with 12 d < P < 40 d and increasing pulsation pe- 
riods as a function of normalized rate of period 
change (upper section). The middle section dis- 
plays running five-point means for the data, and 
the lower section adjusted and extended means of 
the data. 
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Fig. 6. — Normalized blue amplitudes of Cepheids 
with 4 d < P < 8 d and increasing pulsation periods 
as a function of normalized rate of period change 
(upper section). The middle section displays run- 
ning five-point means for the data, and the lower 
section adjusted and extended means of the data 
after judicious separation of the overlapping sam- 
ples using the results of Fig. 5 as a template. 
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Fig. 7. — Normalized blue amplitudes of Cepheids 
with 4 d < P < 8 d and decreasing pulsation pe- 
riods as a function of normalized rate of period 
change (upper section). The middle section dis- 
plays running five-point means for the data, and 
the lower section adjusted and extended means of 
the data after judicious separation of the overlap- 
ping samples as in Fig. 6. 
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